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Abstract — Linear active-power-only DC power flow approx- 
imations are pervasive in the planning and control of power 
systems. However, these approximations fail to capture reactive 
power and voltage magnitudes, both of which are necessary in 
many applications to ensure voltage stability and AC power 
flow feasibility. This paper proposes linear-programming models 
(the LPAC models) that incorporate reactive power and voltage 
magnitudes in a linear power flow approximation. The LPAC 
models are built on a convex approximation of the cosine terms 
in the AC equations, as well as Taylor approximations of the 
remaining nonlinear terms. Experimental comparisons with AC 
solutions on a variety of standard IEEE and MATPOWER 
benchmarks show that the LPAC models produce accurate 
values for active and reactive power, phase angles, and voltage 
magnitudes. The potential benefits of the LPAC models are 
illustrated on two "proof-of-concept" studies in power restoration 
and capacitor placement. 

Index Terms — DC power flow, AC power flow, LP power flow, 
linear relaxation, power system analysis, capacitor placement, 
power system restoration 



Nomenclature 

7 AC Current 

V = v + id AC voltage 
S = p + iq AC power 

Z = r + ix Line impedance 

V = g + ib Line admittance 
Y b = g y + ib y Y-Bus element 
Y c = g c + ib c Line charge 

Y s = g s + ib s Bus shunt 

T = t + is Transformer parameters 

V = \V\Z9° Polar form 

S n AC Power at bus n 

S nm AC Power on a line from n to m 

VM Power network 

N Set of buses in a power network 

L Set of lines in a power network 

G Set of voltage controlled buses 
Slack Bus 

\V h \ Hot-Start voltage magnitude 

| V* | Target voltage magnitude 

cj) Voltage magnitude change 

A Absolute difference 

5 Percent difference 

x Approximation of x 

x Upper bound of x 

x Lower bound of x 
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I. Introduction 

OPTIMIZATION technology is widely used in modern 
power systems JT] and has resulted in millions of dollars 
in savings annually [2]. But the increasing role of demand 
response, the integration of renewable sources of energy 
and the desire for more automation in fault detection and 
recovery pose new challenges for the planning and control of 
electrical power systems |3|. Power grids now need to operate 
in more stochastic environments and under varying operating 
conditions, while still ensuring system reliability and security. 

Optimization of power systems encompasses a broad spec- 
trum of problem domains, including optimal power flow 0J, 

0, i), 0, CD, ED, 03, CH, 02, 03. 03), uvrp-base 

market calculations fl5l . fl6l . ifTTI . transmission switching 
ifTSIl . 03, ED, real-time security-constrained dispatch [21 1, 
ll22ll . day-ahead security-constrained unit commitment [23 1, 
E3), ll25ll . distribution network configuration [26|, [27], capac- 
itor placement [28 1, l|29ll . OUll . expansion planning [31 1, |32|, 
1331 . 1341 . E3, E§1, E3, ED, EH, vulnerability analysis 
BOl . ATI . B2l . 031 . 03), and power system restoration 03), 
BSll to name a few. Some of these use active power only 
while others consider both active and reactive power. 

Restricting attention to active power is often appealing 
computationally as the nonlinear AC power flow equations 
can then be approximated by a set of linear equations that 
define the so-called Linearized DC (LDC) model. Under 
normal operating conditions and with some adjustment for 
line losses, the LDC model produces a reasonably accurate 
approximation of the AC power flow equations for active 
power 03. Moreover, the LDC model can be embedded in 
Mixed-Integer Programming (MIP) models for a variety of 
optimization applications in power system operations. This is 
particularly attractive as the computational efficiency of Linear 
Programming (LP) and MIP solvers has significantly improved 
over the last two decades [48]. 

However, the LDC model does not capture reactive power 
and hence cannot be used for applications such as capacitor 
placement and voltage stability to name only two. Moreover, 
the accuracy of the LDC model outside normal operating 
conditions is an open point of discussion (e.g., |fT5l . 03, |49|, 
ED, BT1 ). This in turn raises concerns for other applications 
such transmission planning, vulnerability analysis, and power 
restoration, which may return infeasible or suboptimal solu- 
tions when the LDC model is used to approximate the AC 
power flow equations. As a result, these applications often 
turn to nonlinear programming techniques (8), 03, 03), 
03, iterative heuristics and decomposition IfTTI . Ifl2l . [36|, 
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[42 1, model relaxation Q, 11281 . tabu search 11301 . and genetic 
algorithms 11291 . ll34ll to ensure feasibility. These techniques 
often require extensive tuning for each problem domain, 
may consume significant computational resources, and cannot 
guarantee global optimality. 

This paper aims at bridging the gap between the LDC model 
and the AC power flow equations. It presents linear programs 
to approximate the AC power flow equations. These linear 
programs, called the LPAC models, are based on two ideas: 

1) They reason both on the voltage phase angles and the 
voltage magnitudes, which are coupled through equa- 
tions for active and reactive power; 

2) They use a piecewise linear approximation of the cosine 
term in the power flow equations and Taylor series for 
approximating the remaining nonlinear terms. 

The LPAC models have been evaluated experimentally over 
a number of standard benchmarks under normal operating 
conditions and various contingencies. Experimental compar- 
isons with AC solutions on standard IEEE and MatPower 
benchmarks shows that the LPAC models are highly accurate 
for active and reactive power, phase angles, and voltage mag- 
nitudes. Moreover, the LPAC models can be integrated in MIP 
models for applications reasoning about reactive power (e.g., 
capacitor placement) or topological changes (e.g., transmission 
planning, vulnerability analysis, and power restoration). 

This rest of this paper presents a rigorous and systematic 
derivation of the LPAC models, experimental results about 
their accuracy, and its application to power restoration and 
capacitor placement. Section [II] reviews the AC power flow 
equations. Section III derives the LPAC models and Section 



form, the real (p n ) and imaginary (q n ) terms become 



|rV] presents the experimental results on its accuracy. Sec- 
tion [V] presents the "proof -of-concept" experiments in power 
restoration and capacitor placement to demonstrate potential 



applications of the LPAC models. Section VI discusses related 



^ \V„\ 2 g nm - \V n \\Vm\{g nm cos{6° n -6° m )+b nm sin(0°-0£J) 

J 7 i 

22 -\V„\ 2 b nm - \V n \\V m \(g nm s\n(8° n -Q° m )-b nrn cos(0°— 6° m )) 

m 

The Y-Bus Matrix: The formulation can be simplified further 
by using a Y-Bus Matrix, i.e., a precomputed lookup table 
that allows the power flow at each bus to be written as a 
summation of 2n terms instead a summation of 3(n— 1) terms. 
Observe that, in the power flow equations ([lj, the first term 
V n V*Y* is a special case of the second term —V n Vj^Y* m 
with — V m — V n . We can eliminate this special case by 

1) extending the summation to include n terms, i.e., ^2^f m 
becomes ^2,', 

2) defining the Y-Bus admittance Y^ m as 



Y" 



Y° — -Y 

1 nm 1 n 



Given the Y-Bus, the power flow equations ([T]i can be rewritten 
as a single summation 



work and Section VII concludes the paper. 



m 

giving us the popular formulation of active and reactive power: 

Pn =Y)V n \\V m M m cos(C-C)+^ m sin(0°-O (3) 

m 

In =Y)Vn\\V m \{9 V nrn Sm(9° n - d° m ) - COs(0°- O (4) 



II. Review of AC Power Flow 
The steady state AC power for bus n is given by 

°n — / t 'n'n^nm y ri v m 1 nm- 



(1) 



This equation is not symmetric. From the perspective of bus 
n, the power flow on a line to bus n is 



nm 



y y*y* _yy*y 

y n * n 1 nm v n v rn 1 1 

while, from the perspective of bus m, it is 



Y V Y* — V V Y 

Vmv m 1 mn v my n 1 mn- 



In general, S nm ^ S mn . 



A. The Traditional Representation 

The AC power flow definition is typically expanded in terms 
of real numbers only. By representing power in rectangular 



B. An Alternate Representation 

The Y-Bus formulation is concise but makes it difficult to 
reason about the power flow equations. This paper uses the 
more explicit equations which can be presented as bus and 
line equations as follows: 



Pn 

q n 

Prim 
Qnm 



^ ^ Prim, 
m 

^ Qnm 

m 

- \V n \ 2 g nm - \V n \\V m \g nm cos(6>° - 6° n ) 

-\v n \\v m \b nm sm{ei-e° m ) 

■ -\V n \\ im + \V n \\V m \b nm cos(0° - C) 
-iKHKilflWi sin(6>° - 9° m ) 



(5) 
(6) 

(7) 
(8) 



Once again, [7] and [8] are asymmetric and the line admittance 
values Y have not been modified. 
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C. Extensions for Practical Power Networks 

In the above derivation, each line is a conductor with an 
impedance Z. The formulation can be extended to line charg- 
ing and other components such as transformers and bus shunts, 
which are present in nearly all AC system benchmarks. We 
show how to model these extensions in the Y-Bus formulation 
for simplicity. 

Line Charging: A line connecting buses n and m may 
have a predefined line charge Y c . Steady state AC models 
typically assume that a line charge is evenly distributed across 
the line and hence it is resonable to assign equal portions of 
its charge to both sides of the line. This is incorporated in the 
Y-Bus matrix as follows: 

yb' yb i yc /o 
nn nn 1 rival > 

yb' = yb i y-c /o 
mm mm 1 nm/ 

Transformers: A transformer connecting bus n to bus 
m can be modeled as a line with modifications to the Y-Bus 
matrix. The properties of the transformer are captured by a 
complex number T nm = \T\ Z.s° , where |T| is the tap ratio 
from n to m and s° is the phase shift. It is worth noting 
that the direction of a transformer-line is very important to 
model the tap ratio and phase shift properly. A transformer is 
modeled in the Y-Bus matrix as follows: 

yb' —Y b —Y -\-Y /IT I 2 
1 nn 1 nn 1 nm \ 1 nml\- L nm\ ) 

yb' y Irp* 

1 nm ~ 1 nm/ J-nnn 

yb' _ y Iji 
mn — mn j - 1 - nm ■ 

If a line charge exists, it must be applied before the transformer 
calculation, i.e., 

Y nn — Y nn — Y nm + (Y nm + Y nm / 2) / \T nm \ 

Bus Shunts: A bus n may have a shunt element which is 
modeled as a fixed admittance to ground with a value of Y s . 
In the Y-Bus matrix, we have 

yb' _ yb , ys 

nn nn ' n 

Unlike line charging, this extension is not affected by trans- 
formers, since it applies to a bus and not a line. 

D. The Linearized DC Power Flow 

Many variants of the Linearized DC (LDC) model exist 
1 52 1, 1531 , 1541 , 1531 . A comprehensive review of all these 
variants is outside the scope of this work but an in-depth 
discussion can be found in [47). For brevity, we only review 
the simplest and most popular variant of the LDC, which is 
derived from the AC equations through a series of approxi- 
mations justified by operational considerations under normal 
operating conditions. In particular, the LDC assumes that (1) 
the susceptance is large relative to the conductance |<7| <C \b\; 
(2) the phase angle difference is small enough to ensure 
sin(0° - 6>° n ) « 9° — 9 m ; and (3) the voltage magnitudes 
| V| are close to 1.0 and do not vary significantly. Under these 
assumptions, Equations (|7]i and |8]l reduce to 

Pnm = -bnm(0° n -6 m ) (9) 



This simple linear formulation has been used in many frame- 
works for decision support in power systems lfT31 . lfT8l . OTl . 
EH, P31 . Il46l . This traditional model is used as the baseline 
in the experimental results. 

III. Linear-Programming Approximations 

This section presents linear-programming approximations 
of the AC power flow equations. To understand the approxi- 
mations, it is important to distinguish between hot-start and 
cold-start contexts (47] ■ In hot-start contexts, a solved AC 
base-point solution is available and hence the model has at its 
disposal additional information such as voltage magnitudes. 
In cold start contexts, no such solved AC base-point solution 
is available and it can be "maddeningly difficult" |[T31 to 
obtain one by simulation of the network. Hot-start models are 
well-suited for applications in which the network topology is 
relatively stable, e.g., in LMP-base market calculations, op- 
timal line switching, distribution configuration, and real-time 
security constrained economic dispatch. Cold-start models are 
used when no operational network is available, e.g., in long- 
term planning studies. We also introduce the concept of warm- 
start contexts, in which the model has at its disposal target 
voltages (e.g., from normal operating conditions) but an actual 
solution may not exist for these targets. Warm-start models 
are particularly useful for power restoration applications in 
which the goal is to return to normal operating conditions as 
quickly as possible. This section presents the hot-start, warm- 
start, and cold-start models in stepwise refinements. It also 
discusses how models can be generalized to include generation 
and load shedding, remove the slack bus, impose constraints 
on voltages and reactive power, and capacity constraints on 
the lines, all which are fundamental for many applications. 

A. AC Power Flow Behavior 

Before presenting the models, it is useful to review the 
behavior of AC power flows, which is the main driver in 
the derivation. The high-level behavior of power systems is 
often characterized by two rules of thumb in the literature: (1) 
phase angles are the primary factor in determining the flow 
active power; (2) differences in voltage magnitudes are the 
primary factor in determining the flow of reactive power [56|. 
We examine these properties experimentally. 

The experiments make two basic assumptions: (1) In the 
per unit system, voltages do not vary far from a magnitude of 
1.0 and angle of 0.0; (2) The magnitude of a line conductance 
is much smaller than the magnitude of the susceptance, i.e., 
\g\ <C |6|. We can then explore the bounds of the power flow 
equations Jfy and ([HJl^when the voltages are in the following 
bounds: \V n \ = l.0,\V m \ G (1.2, O.8),0°-0^ € (-tt/6, tt/6). 
These bounds are intentionally generous so that the power flow 
behavior within and outside normal operating conditions may 
be illustrated. 

Figure [T] presents the contour of the active power (left) 
and reactive power (right) equations for a line (n, m) under 
these assumptions when Y nm — 0.2 — tl. The contour lines 
indicate significant changes in power flow. Consider first the 
active power plot (left). For a fixed voltage, varying the phase 
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Active Power Field 



Reactive Power Field 





Voltage Difference 



Fig. 1. Power Flow Contour of Active (left) and Reactive (right) Power with 
g = 0.2 and b = -1. 

angle difference induces significant changes in active power 
as many lines are crossed. In contrast, for a fixed phase angle 
difference, varying the voltage has limited impact on the active 
power, since few lines are crossed. Hence, the plot indicates 
that phase angle differences are the primary factor of active 
power flow while voltage differences have only a small effect. 
The situation is quite different for reactive power (right plot). 
For a fixed voltage, varying the phase angle difference induces 
some significant changes in reactive power as around four 
lines can be crossed. But, if the phase angle difference is 
fixed, varying voltage induces even more significant changes 
in reactive power since as many as seven lines may now be 
crossed. Hence, changes in voltages are the primary factor of 
reactive power flows but the phase angle differences also have 
a significant influence. 

B. The Hot-Start LPAC Model 

The linear-programming approximation of the AC Power 
flow equations in a hot-start context is based on three ideas: 

1) It uses the voltage magnitude \V^\ from the AC base- 
point solution at bus n; 

2) It approximates sin(a;) by x; 

3) It uses a convex approximation of the cosineQ 

Let cos(0° — f?£j denote the convexification of cos(f?° — 0£J in 
the range (— 7r/2,7r/2), then the linear-programming approx- 
imation solves the line flow constraints 



i I llll 



g nm cos( 



\v£\\v*\b 

\V£\ 2 b nm 

"rht iTrh 



W„ h \\V*\b n 



V:\\V»\g nm (6° n -e°J 



(10) 



(11) 



The details of the convexification are given in Appendix [A] 
The hot-start model, presented in Model [T] is a linear program 
that replaces the AC power equations by Equations 10 and 11 
and thus captures an approximation of reactive power in a 
linear formulation. 

A complete linear program for this formulation is presented 
in Model [T] The inputs to the model are: (1) A power network 
VM = (N, L, G, s), where N is the set of buses, L is the set 
of lines, G is the set of voltage-controlled generators, s is 

'The domain of the cosine should not exceed the range (— 7r/2,7r/2) to 
ensure convexity. This range is generous for AC power flows. 



Model 1 The Hot-Start LPAC Model. 
Inputs: 

VM = (N, L, G, s) 
\V h \ 



phase angle on bus n (radians) 
Approximation of cos(0° — 6*^ ) 



Vn G N n ^ s 
VneJV n ^ s ri<^ G 



(N 1 

(MI 

(IV 1 

(N 1 



the power network 

voltage magnitudes from a base-point solution 
cs - cosine approximation segment count 

Variables: 

8° e (-00,00) 

COSnm 6 (0, 1) 

Maximize: 

(n,m) £2. 

gSubject to: 

e° = o 

mSiV 

4n = «L 

mSiV 

V(n, m), (m, n) 6 L 

: \Vn\ 2 9nrn-\V£\\V£\(g nm COS nm + b nm {9 



v h = 



\v»\ 2 b nm -\v£\\vMg nm (e° - e° m ) - b nm £Ss nm ) (tvi 



PWL(COS}(cos nm ,(6>° 



J,— 7r/3,7r/3, cs) 



5) 
6) 

(Mr 7) 



the slack bus; (2) the voltage magnitudes \V \ for the buses 
and (3) the number of segments cs for approximating the 
cosine function. The objective (fV(T] 1 ) maximizes the cosine 
approximation to make it as close as possible to the true 
cosine value. Constraints (IvJT]2) model the slack bus, which 
has a fixed phase angle. Constraints (1V(T]3) and (M(T|4) model 
KCL on the buses. Like in AC power flow models, the KCL 
constraints are not enforced on the slack bus for both active 
and reactive power and on voltage-controlled generators for 
reactive power. Constraints (N^T]5) and (N^T|6) capture the 
approximate line flows from Equations ( JT0) > and (pT[ ). Finally, 
Constraints (1V|T|7) define a system of inequalities capturing 
the piecewise-linear approximation of the cosine terms in the 
domain (— 7r/3,7r/3) using cs line segments for each line in 
the power network. 

To our knowledge, this hot-start model is the first linear 
formulation that captures the cosine contribution to reactive 
power. However, fixing the voltage magnitudes, \V h \, in the 
power flow equations may be too restrictive in many applica- 
tions. In the remaining sections, we remove this restriction. 

C. The Warm-Start LPAC Model 

This section derives the warm-start LPAC model, i.e., a 
Linear-Programming model of the AC power flow equations 
for the warm-start context. The warm-start context assumes 
that some target voltages | V"* | are available for all buses except 
voltage-controlled generators whose voltage magnitudes \V 9 \ 
are known. The network must operate close (e.g., ±0.1 Volts 
p.u.) to these target voltages, since otherwise the hardware 
may be damaged or voltages may collapse. 

The warm-start LPAC model is based on two key ideas: 

1) The active power approximation is the same as in the 
hot-start model, with the target voltages replacing the 
voltages in the base-point solution; 

2) The reactive power approximation reasons about voltage 
magnitudes, since changes in voltages are the primary 
factor of reactive power flows. 
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To derive the reactive power approximation in the warm-start 
LPAC model, let <fi be the difference between the target voltage 
and the true value, i.e., 

\v\ = \v t \ + ct>. 

Substituting in Equation [8] we obtain 

{\KW m \ + \V*\<f>m + \v^\4> n + 4>M 

(g nm sin(6° n - e°J - b nm cos(6 n - o (12) 

We can divide this expression into two parts 



Qnm — Qnm + Inm 



(13) 



where q t nm is Equation [i] with \V\ = \V t \ and q^ m captures 



the remaining terms, i.e., 



-(2\v*\<p n +yl)b nm - 
(I^|0 m + |v;*j0 n + ^ m ) 

(g nm sin(6»° - 9° m ) - b nm cos(fr'„ 



(0°-O (14) 



Equation 13 is equivalent to Equation [8] and must be linearized 
to obtain the LPAC model. 

The q t nm part has target voltages and may thus be approx- 
imated like q\ m . The q^ m is more challenging as it contains 
nonlinear and non-convex terms such as <fi n 4>m cos(6>° — f?£j. 
We approximate q£ m using the linear terms of the Taylor series 
of q£ m at (j) n = 0, 4> m = 0, 6° — 8° n = to obtain 



q* m = -(2\Vt\<t> n )b„ 



(\V^ rn + \Vl\(b n )b nm (15) 



or, equivalently, 



turn = -\V*Km{ct> n ~ 4> m ) ~ {K\ ~ Wn (16) 

A complete linear program for this formulation is presented in 
Model [2] The inputs^ to the model are similar to ModejJT] with 
hot start voltages \V h \ replaced by target voltages \V t \. The 
objective (N^2]l) maximizes the cosine approximation to make 
it as close as possible to the true cosine value. Constraints 
(]Vj2]2) model the slack bus, which has a fixed voltage and 
phase angle. Constraints (1V(2]3) capture the voltage-controlled 
generators which, by definition, do not vary from their voltage 
target Constraints (M2]4) and (M5]5) model KCL on 

the buses, as well as the effects of voltage change presented 
in Equation ( [16) . Like in AC power flow models, the KCL 
constraints are not enforced on the slack bus for both active 
and reactive power and on voltage-controlled generators for 
reactive power. Constraints (IV|2]6) and (M5]7) capture the 
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approximate line flows from Equations 
straints (N^SjS) model the effects of voltage change presented 
in Equation ( p"6] ). Finally, Constraints (N^2]9) define a system 
of inequalities capturing the piecewise-linear approximation 
of the cosine terms in the domain (— 7r/3,7r/3) using cs line 
segments for each line in the power network. 



Model 2 The Warm-Start LPAC Model. 
Inputs: 

VN - 



\v*\ 

cs 



(N,L,G,s) 



Variables: 

o° e (-00,00) 

e (— |V*|,O0) 

COSnm £ (0, 1) 

Maximize: 



the power network 
target voltage magnitudes 
cosine approximation segment count 

phase angle on bus n (radians) 
voltage change on bus n (Volts p.u.) 
Approximation of cos(0° — 8^) 



E 



Vn G N n ^ s 



(n,m) £ L 

Subject to: 

6° = 0,(j> s = 

4>i = Vi e G 

Vn = Pnm 

mSiV 

In = y t Qnm 

V(n, m), (m, n) £ L 

■-\yi?9nm-\ K\ I % I ( 9 „ m cos nm + b nm (0\ 

-\vl\ 2 b nm -\v^\\v^\{g nrn {e° n - e° m ) - b. 

= -\V^\b nm {K - 4>m) - (|V^| - \V m \)b 
PWL(C05>(cos nm , (6° -9° m ), -tt/3, tt/3, cs) 



(V2 



(IV 2 



2) 
(M2 3) 



inm 



Vn 6 N n ^ s n £ G 



(V2 



(IV 2 



4) 



5) 



Pnm 
<7* = 

■Inm 



T,COS nm ) 



(V2 6) 
(V2 7) 
(IV 2 8) 
(IV 2 9) 



D. The Cold-Start LPAC Model 

We now conclude by presenting the cold-start LPAC model. 
In a cold-start context, no target voltages are available and 
voltage magnitudes are approximated by 1.0, except for 
voltage-controlled generators whose voltages are given by \V%\ 
(n E G). The cold-start LPAC model is then derived from the 
warm-start LPAC model by fixing \V*\ = 1 for all i G N. 



Equation (|T6jl then reduces to 

= -b n , 



q A 



1) 



(17) 



Figure [3] presents the cold-start LPAC model, which is very 
close to the warm-start model. Note that Constraints (1V(3]3) 
use 4>i to fix the voltage magnitudes of generators. 

E. Extensions to the LPAC Model 

The LPAC models can be used to solve the AC power flow 
equations approximately in a variety of contexts. This section 
reviews how to generalize the LPAC models for applications 
in disaster management, reactive voltage support, transmission 
planning, and vulnerability analysis. The extensions are illus- 
trated on the warm-start model but can be similarly applied to 
the cold-start model. 

Generators: The LPAC model can easily be generalized 
to include ranges for generators: Simply remove the generator 
from G and place operating limits on the p and q variables 
for that bus. In this formulation, voltage-controlled generators 
can also be accommodated by fixing <p n to zero at bus n. 

Removing the Slack Bus: By necessity, AC solvers use a 
slack bus to ensure the flow balance in the network when the 
total power consumption is not known a priori (e.g., due to line 
losses). As a consequence, the LPAC model depicted in Figure 
[2] also uses a slack bus so that the AC and LPAC models can 
be accurately compared in our experimental results. However, 
it is important to emphasize that the LPAC model does not 
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Model 3 The Cold-Start LPAC Model 
Inputs: 

VM = (TV, L, G, s) 
cs 

Variables: 

9° G (-00,00) 



the power network 
cosine approximation segment count 



phase angle on bus ri (radians) 
voltage change on bus n (Volts p.u.) 
Approximation of cos(0° — 9 m ) 



_e(-|v*|,oo) 

COSnm £ (0, 1) 

Maximize: 

^ ^ COSnm 
(n,m) £L 

Subject to: 

9° = 0, = 

<t>i = \K 3 \ ~ L0 Vi e G 

Pn= ^2 P nm Vn 6 N 11 ^ S 

9™ = 5nm + 9nm Vn 6 JV n^s n^G 

mSiV 

V(n, m), (m, n) £ L 

~ 9nm{9 n 9 m ) b nm cos nrn 



(TV 3 



(TV 3 
(V3 

(TV 3 



(TV 3 



5) 



1 ?"i /?), 

(5* = 

<?nm = -Km{<j>n 

PWL(C05>(cos I , 



x) 



,),-7r/3,7r/3,cs) 



(IV 3 6) 
(VI 7) 

(V58) 
(V3"9) 



need a slack bus and the only reason to include a slack bus in 
this model is to allow for meaningful comparisons between the 
LPAC and AC models. As discussed above, the LPAC model 
can easily include a range for each generator, thus removing 
the need for a slack bus. 

Load Shedding: For applications in power restoration 
(e.g., 03), 06), ED), the LPAC model can also integrate load 
shedding: Simply transform the loads into decision variables 
with an upper bound and maximize the load served. The cosine 
maximization should also be included in the objective but with 



a smaller weight. Section V-A reports experimental results on 



such a power restoration model. 

Modeling Additional Constraints: In practice, feasibility 
constraints may exist on the acceptable voltage range, the 
reactive injection of a generator, or line flow capacities. 
Because Model [2] is a linear program, it can incorporate such 
constraints. For instance, constraint 

\V\<\V^\+4> n VneiV 

ensures that voltages are above a certain limit \V\, constraint 

limits the maximum reactive injection bounds at bus n to q^. 
Finally, let \S nm \ be the maximum apparent power on a line 
from bus n to bus m. Then, constraint 



(Pnm) + i'lnm + 9nm) <l^nro| 

ensures that line flows are feasible in the LPAC model. The 
quadratic functions can be approximated by piecewise-linear 
constrains (e.g., Il50l ). 



IV. Accuracy of the LPAC Model 

This section evaluates the accuracy of the LPAC models 
by comparing them to an ideal nonlinear AC power flowj^] 
It includes a detailed analysis of the model accuracy (Sec- 



tion [TV-Ai and an investigation of alternative approximations 
(Section IV-B i. The experiments were performed on nine 
traditional power-system benchmarks which come from the 
IEEE test systems E) and Matpower J58). The AC power 
flow equations were solved with a Newton-Raphson solver 
which was validated using MATPOWER. The LPAC models 
use 20 line segments in the cosine approximation and all 
of the models solved in less than 1 second on a 2.5 GHz 
Intel processor. The results also include a modified version of 
the IEEEddl7 benchmark, called IEEEddl7m. The original 
IEEEddl7 has the slack bus connected to the network by 
a transformer with |T| = 1.05. The nonlinear behavior of 
transformers induces some loss of accuracy in the LPAC model 
and, because this error occurs at the slack bus in IEEEddl7, it 
affects all buses in the network. IEEEddl7m resolves this issue 
by setting |T| = 1.00 and the slack bus voltage to 1.05. As 
the results indicate, this equivalent formulation is significantly 
better for the LPAC model. 

A. Accuracy of The LPAC Models 

This section reports empirical evaluations of the LDC and 
LPAC models in cold-start and warm-start contexts. It reports 
aggregate statistics for active power (Table [I]), bus phase angles 
(Table ID), reactive power (Table III i, and voltage magnitudes 
(Table |IV]>. Data for the LDC model is necessarily omitted 



from Tables HI and IV as reactive power and voltages are not 
captured by that model. In each table, two aggregate values 
are presented: Correlation (corr) and absolute error (A). The 
units of the absolute error are presented in the headings. Both 
average (//) and worst-case (max) values are presented. The 
worst case can often be misleading: For example a very large 
value may actually be a very small relative quantity. For this 
reason, the tables show the relative error (<5) of the value 
selected by the max operator using the arg-max operator. The 
relative error is a percentage and is unit-less. 

Table [I] indicates uniform improvements in active power 
flows, especially in the largest benchmarks IEEE118, 
IEEEddl7, and MP300. Significant errors are not uncommon 
for the linearized DC model on large benchmarks ll47l and are 
primarily caused by a lack of line losses. Due to its asymmet- 
rical power flow equations and the cosine approximation, the 
LPAC model captures line losses. 

Table[n|presents the aggregate statistics on bus phase angles. 
These results show significant improvements in accuracy es- 
pecially on larger benchmarks. The correlations are somewhat 
lower than active power, but phase angles are quite challenging 
from a numerical accuracy standpoint. 



Table III presents the aggregate statistics on line reactive 
power flows. They indicate that reactive power flows are gen- 
erally accurate and highly precise in warm-start contexts. To 
highlight the model accuracy in cold-start contexts, the reactive 

2 For consistency, the LPAC models are extended to include line charging, 
bus shunts, and transformers, as discussed in Section III-CI 
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TABLE I 

Accuracy of the LPAC Model: Active Power Flows. 



TABLE II 

Accuracy of the LPAC Model: Phase Angles. 



Benchmark 




Active Power (MW) 




Corr 


M(A) 


max(A) 


(5(arg-max(A)) 


The LDC Model 


ieeel4 


0.9994 


1.392 


10.64 


6.783 


mp24 


0.9989 


5.659 


19.7 


23.65 


ieee30 


0.9993 


1.046 


13.1 


7.562 


mp30 


0.9993 


0.2964 


2.108 


19.36 


mp39 


0.9995 


7.341 


43.64 


6.527 


ieee57 


0.9989 


1.494 


8.216 


8.055 


ieeel 18 


0.9963 


3.984 


56.3 


44.74 


ieeeddl 7 


0.9972 


4.933 


201.3 


13.84 


ieeeddl 7m 


0.9975 


4.779 


191.1 


13.23 


mp300 


0.9910 


11.09 


418.5 


90.02 


The LPAC-Cold Model 


ieeel 4 


0.9989 


1.636 


5.787 


13.13 


mp24 


0.9999 


1.884 


6.159 


2.933 


ieee30 


0.9998 


0.5475 


2.213 


2.523 


mp30 


0.9995 


0.2396 


1.641 


15.07 


mp39 


1.0000 


2.142 


8.043 


3.288 


ieee57 


0.9995 


0.9235 


4.674 


9.728 


ieee 118 


1.0000 


0.622 


3.708 


2.038 


ieeeddl 7 


0.9999 


1.827 


30.38 


2.088 


ieeeddl 7m 


0.9999 


1.475 


20.21 


1.399 


mp300 


0.9998 


2.455 


18 


8.675 


The LPAC-Warm Model 


ieee 14 


1.0000 


0.1689 


1.588 


1.012 


mp24 


1.0000 


0.6621 


2.041 


1.01 


ieee30 


1.0000 


0.1847 


2.433 


1.405 


mp30 


0.9999 


0.1052 


0.705 


6.474 


mp39 


1.0000 


1.557 


11.58 


1.731 


ieee57 


1.0000 


0.2229 


2.013 


1.973 


ieee 118 


0.9999 


0.4386 


7.376 


5.862 


ieeeddl 7 


1.0000 


0.58 


22.5 


1.547 


ieeeddl 7m 


1.0000 


0.5725 


21.73 


1.504 


mp300 


0.9999 


1.195 


52.84 


11.37 



Line Reactive Power Correlation (MVar) 



Line Reactive Power Correlation (MVar) 




AC Power Flow 



Benchmark 




Phase Angle (rad) 




Corr 


MA) 


max(A) 


<5(arg-max(A)) 


The LDC Model 


ieee 14 


0.9993 


0.02487 


0.04258 


15.22 


mp24 


0.9997 


0.01334 


0.02037 


15.23 


ieee30 


0.9981 


0.02831 


0.04733 


16.45 


mp3() 


0.9800 


0.005658 


0.01607 


30.27 


mp39 


0.9951 


0.0283 


0.05813 


85.56 


ieee57 


0.9898 


0.02244 


0.05958 


24.1 


ieeel 18 


0.9904 


0.03452 


0.09026 


88.41 


ieeeddl 7 


0.9892 


0.115 


0.1395 


16.09 


ieeeddl 7m 


0.9920 


0.0461 


0.06924 


41.88 


mp300 


0.9752 


0.3103 


0.4244 


975.7 


The LPAC-Cold Model 


ieee 14 


0.9971 


0.004525 


0.01241 


5 


mp24 


0.9999 


0.003539 


0.008947 


6.922 


ieee30 


0.9965 


0.007268 


0.02413 


8.386 


mp30 


0.9782 


0.006236 


0.01804 


33.99 


mp39 


0.9989 


0.006268 


0.02314 


34.06 


ieee57 


0.9894 


0.0179 


0.05467 


22.11 


ieee 118 


0.9994 


0.003225 


0.01354 


9.633 


ieeeddl 7 


0.9981 


0.03648 


0.05165 


5.958 


ieeeddl 7m 


0.9990 


0.007207 


0.02682 


4.522 


mp300 


0.9984 


0.01458 


0.08086 


38.49 


The LPAC-Warm Model 


ieee 14 


1.0000 


0.001448 


0.001829 


0.6914 


mp24 


1.0000 


0.001337 


0.002203 


2.156 


ieee30 


1.0000 


0.002345 


0.002819 


0.9629 


mp3() 


0.9998 


0.001298 


0.001774 


4.775 


mp39 


0.9999 


0.005315 


0.006241 


4.273 


ieee57 


1.0000 


0.002711 


0.00357 


1.776 


ieee 118 


0.9999 


0.005958 


0.008366 


2.526 


ieeeddl 7 


0.9999 


0.01492 


0.01719 


2.419 


ieeeddl 7m 


0.9999 


0.008443 


0.01059 


2.33 


mp300 


0.9997 


0.03842 


0.04502 


18.51 



Bus Voltage Correlation (Volts p.u.) 






Bus Voltage Correlation {Volts p 


u.) 



















































AC Powar Flow 



Fig. 2. Reactive Power Flow Correlation for the LPAC Model on IEEEddl7m Fig. 3. Voltage magnitude correlation for Model LPAC on IEEEddl7m (left) 
(left) and MP300 (right) in a cold-start context. and MP300 (right) in a cold-start context. 



flow correlation for the two worst benchmarks, IEEEddl7m 
and MP300, is presented in Figure [2] 



Table IV presents the aggregate statistics on bus voltage 
magnitudes. These results indicate that voltage magnitudes are 
very accurate on small benchmarks, but the accuracy reduces 
with the size of the network. The warm-start context brings a 
significant increase in accuracy in larger benchmarks. To illus- 
trate the quality of these solutions in cold-start contexts, the 
voltage magnitude correlation for the two worst benchmarks, 
i.e., IEEEddl7m and MP300, is presented in Figure [3] The 
increase in voltage errors is related to the distance from a load 
point to the nearest generator. The linearized voltage model 
incurs some small error on each line. As the voltage changes 



over many lines, these small errors accumulate. By comparing 
the percentage of voltage-controlled generator buses in each 
benchmark |G|/|7Y| (Table |v| to accuracy in Table IV the 
IEEE57 and IEEEddl7 benchmarks indicate that that a low 
percentage is a reasonable indicator of the voltage accuracy in 
the cold-start context. 



B. Alternative Linear Models 

The formulation of the LPAC models explicitly removes two 
core assumptions of the traditional LDC model: 

1) Although cos(6>° — 9° m ) maybe very close to 1, those 
small deviations are important. 
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TABLE III 

Accuracy of the LPAC Model: Reactive Power Flows. 



TABLE V 

Percentage of Voltage-Controlled Buses in the Benchmarks. 



Benchmark 




Reactive Power (MVar) 




Corr 


M(A) 


max(A) 


(5(arg-max(A)) 


The LPAC-Cold Model 


ieeel4 


0.9948 


0.7459 


2.561 


14.92 


mp24 


0.9992 


1.505 


5.245 


9.309 


ieee30 


0.997 


0.4962 


1.902 


23.36 


mp30 


0.9991 


0.3135 


0.8925 


3.886 


mp39 


0.9973 


3.898 


15.15 


18.25 


ieee57 


0.9991 


0.5316 


2.98 


3.973 


ieeel 18 


0.9991 


0.7676 


6.248 


8.561 


ieeeddl7 


0.9789 


3.989 


48.65 


40.58 


ieeeddl7m 


0.9927 


2.415 


34.18 


12.36 


mp300 


0.9981 


3.85 


62.32 


17.98 


The LPAC-Warm Model 


ieeel 4 


0.9895 


0.8689 


3.167 


43.89 


mp24 


0.9992 


1.505 


5.245 


9.309 


ieee30 


0.9975 


0.3455 


1.607 


7.62 


mp30 


0.9991 


0.3135 


0.8925 


3.886 


mp39 


0.9971 


4.03 


15.75 


18.98 


ieee57 


0.9995 


0.3853 


1.46 


5.67 


ieee 118 


0.9992 


0.6326 


6.109 


6.808 


ieeeddl7 


0.9791 


3.985 


48.37 


40.34 


ieeeddl7m 


0.9927 


2.409 


33.96 


12.28 


mp300 


0.9943 


3.584 


162 


45.05 



TABLE IV 

Accuracy of the LPAC Model: Voltage Magnitudes. 



Benchmark 




Voltage Magnitude (Volts p.u.) 




Corr 


M(A) 


max(A) 


5(arg-max(A)) 


The LPAC-Cold Model 


ieee 14 


0.9828 


0.003524 


0.01304 


1.236 


mp24 


0.9983 


0.000676 


0.003244 


0.3362 


ieee30 


0.9908 


0.002445 


0.01098 


1.098 


mp30 


0.9884 


0.002186 


0.009453 


0.9723 


mp39 


0.9992 


0.0007521 


0.002446 


0.2313 


ieee57 


0.9726 


0.01038 


0.03353 


3.587 


ieeel 18 


0.9989 


0.000717 


0.00476 


0.4926 


ieeeddl7 


0.9651 


0.01376 


0.03345 


3.424 


ieeeddl7m 


0.9815 


0.00647 


0.01579 


1.659 


mp300 


0.9948 


0.002361 


0.01552 


1.656 


The LPAC-Warm Model 


ieee 14 


0.9998 


0.0005479 


0.001173 


0.1111 


mp24 


0.9996 


0.000542 


0.002214 


0.2294 


ieee30 


0.9994 


0.001426 


0.002508 


0.25 


mp30 


1.0000 


0.0003884 


0.000707 


0.073 


mp39 


0.9983 


0.00154 


0.003545 


0.3524 


ieee57 


0.9987 


0.002138 


0.005515 


0.59 


ieee 118 


0.9999 


0.0001961 


0.001303 


0.1344 


ieeeddl7 


0.9858 


0.01204 


0.02597 


2.796 


ieeeddl7m 


0.9760 


0.009819 


0.02037 


2.067 


mp300 


0.9967 


0.002477 


0.01403 


1.52 



ieee 14 


mp24 


ieee30 


mp30 


mp39 


35.7% 


45.8% 


20.0% 


20.0% 


30.8% 


ieee57 


ieeel 18 


ieeeddl7 


mp300 




12.3% 


45.8% 


7.4% 


24.7% 





TABLE VI 

Accuracy Comparison of Various Linear Models (Part I). 



Model 


Cumulative Absolute Error 

9t(y n — V m ) \$Z(V n — V m )\ p n 1 q n 


ieee 14 


LDC 


0.3839 


0.166 


13.39 


118.4 


LPAC-A-GC 


0.1561 


0.1296 


13.39 


140.7 


LPAC-A-G 


0.1221 


0.1229 


13.39 


120.7 


LPAC-A-C 


0.1277 


0.1262 


8.843 


53.24 


LPAC 


0.1008 


0.1234 


1.783 


11.43 


mp24 


LDC 


0.448 


0.1434 


53.22 


792.4 


LPAC-A-GC 


0.3676 


0.1289 


53.22 


546.2 


LPAC-A-G 


0.2309 


0.1332 


53.22 


314.3 


LPAC-A-C 


0.2417 


0.116 


53.12 


476.2 


LPAC 


0.03411 


0.0828 


6.94 


64.12 


ieee30 


LDC 


0.5429 


0.2934 


17.55 


169.9 


LPAC-A-GC 


0.1607 


0.2377 


17.55 


173.8 


LPAC-A-G 


0.1284 


0.2268 


17.55 


147.7 


LPAC-A-C 


0.1587 


0.1638 


15.99 


66.28 


LPAC 


0.1305 


0.1476 


2.9 


16.72 


mp30 


LDC 


0.4341 


0.1728 


2.444 


181 


LPAC-A-GC 


0.166 


0.1584 


2.444 


20.78 


LPAC-A-G 


0.1616 


0.1581 


2.444 


17.96 


LPAC-A-C 


0.06886 


0.1454 


2.444 


9.387 


LPAC 


0.07338 


0.1456 


1.736 


6.76 


mp39 


LDC 


0.5634 


0.1997 


43.64 


2816 


LPAC-A-GC 


0.2449 


0.2105 


43.64 


950.3 


LPAC-A-G 


0.1418 


0.1958 


43.64 


284.1 


LPAC-A-C 


0.1896 


0.1868 


37.65 


947.9 


LPAC 


0.04462 


0.1337 


0.4745 


93.31 



two interesting points. First, all linear models tend to bring 
improvements over a traditional LDC model. Second, although 
integrating either the g value or the cosine term brings some 
small improvement independently, together they make signifi- 
cant improvements in accuracy. Additionally a comparison of 
Table [Vl] and Table IVHI reveals that the benefits of the new 
linear models are more pronounced as the network increases. 



2) Although |<jr| <C |£>|, the conductance contributes signifi- 
cantly to the phase angles and voltage magnitudes. 

This section investigates three variants of the cold-start LPAC 
model that reintegrate some of the assumptions of the LDC 
model. The new models are: (1) the LPAC-C model where only 
the cosine approximation is used and g = 0; (2) the LPAC-G 
model where only the g value is used and cos(x) = 1; (3) 
the LPAC-CG model where cos(x) = 1 and g = 0. Tables 



VI and VII present the cumulative absolute error between the 



proposed linear formulations and the true nonlinear solutions. 
Many metrics may be of interest but these results focus on 
line voltage drop V n — V m and bus power S n . These were 
selected because they are robust to errors which accumulate 
as power flows through the network. The results highlight 



V. Case Studies 

This section describes two case studies to evaluate the 
potential of the LPAC models: Power restoration and capacitor 
placement. The goal is not to present comprehensive solutions 
for these two complex problems, but to provide preliminary 
evidence that the LPAC models may be useful in striking a 
good compromise between efficiency and accuracy for such 
applications. This section should be viewed as presenting a 
"proof-of-concept" that the LPAC models may be valuable 
for certain classes of applications where the LDC model is 
not accurate enough and existing approaches are too time 
consuming or suboptimal. 
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TABLE VII 

Accuracy Comparison of Various Linear Models (Part II). 



Model 




Cumulative Absolute Error 






»(Vn-V TO ) 


S(Vn-V m ) 


Pn 




ieee57 


LDC 


1.343 


0.6803 


27.9 


529.1 


LPAC-A-GC 


0.4158 


0.5489 


27.9 


264 


LPAC-A-G 


0.3238 


0.5311 


27.9 


264.1 


LPAC-A-C 


0.3773 


0.4153 


24.43 


113.8 


LPAC 


0.3647 


0.3854 


4.736 


26.35 


ieee 1 1 8 


LDC 


3.083 


1.239 


132.7 


2152 


LPAC-A-GC 


0.7298 


0.9944 


132.7 


1364 


LPAC-A-G 


0.6502 


0.9929 


132.7 


1194 


LPAC-A-C 


0.4417 


0.8553 


104.3 


750 


LPAC 


0.2625 


0.5252 


0.7279 


142.1 


ieeedd!7 


LDC 


4.144 


3.263 


201.3 


3857 


LPAC-A-GC 


5.783 


4.881 


201.3 


2719 


LPAC-A-G 


4.169 


3.242 


201.3 


616.4 


LPAC-A-C 


4.162 


3.48 


200.5 


2660 


LPAC 


1.135 


1.019 


30.38 


362.1 


ieeedd!7m 


LDC 


3.798 


1.972 


191.1 


3353 


LPAC-A-GC 


5.152 


3.219 


191.1 


2210 


LPAC-A-G 


3.302 


2.111 


191.1 


389.4 


LPAC-A-C 


3.31 


2.118 


190.4 


2146 


LPAC 


0.49 


0.6324 


20.21 


223.7 


mp300 


LDC 


13.76 


4.689 


418.5 


14240 


LPAC-A-GC 


11.2 


5.324 


418.5 


5595 


LPAC-A-G 


9.831 


5.171 


418.5 


1648 


LPAC-A-C 


7.855 


4.403 


348.5 


5434 


LPAC 


0.8699 


1.378 


9.703 


976.4 



A. Power Restoration 

After a significant disruption due to, say, a natural disaster, 
large sections of the power network need to be re-energized. To 
understand the effects of restoration actions, power engineers 
must simulate the network behaviour under various courses of 
action. However, the network is far from its normal operating 
state, which makes it extremely challenging to solve the AC 
power flow equations. In fact, the task of finding an AC 
solution without a reasonable starting point has been regarded 
as "maddeningly difficult" 1151 . The LPAC model studied 
here has the benefit of providing starting values for all the 
variables in the AC power flow problem, unlike the traditional 
LDC which only provides active power values. Furthermore, 
the LPAC model has the additional advantage of supporting 
bounds on reactive generation and voltage magnitudes and 
such constraints are critical for providing feasible solutions 
to the AC power flow. This section illustrates these benefits. 

Before presenting the power-restoration model, it is impor- 
tant to mention the key aspect of this application. When the 
power system undergoes significant damages, load shedding 
must occur. The LDC and LPAC models must be embedded 
in a restoration model that maximizes the served load given 
operational constraints such as the generation limits. These 
load values indicate the maximum amount of power that can be 
dispatched while ensuring system stability. Model [4] presents 
a linear program based on the warm-start LPAC model which, 
given limits on active power generation pS and the desired 
active and reactive loads p l , q l at each bus, determines the 



Model 4 A LP for Maximizing Desired Load. 
Inputs: 

p n - maximum active injection for bus n 

p l n - desired active load at bus n 

q l n - desired reactive load at bus n 

Inputs from Model [2] (The Warm-Start LPAC Model) 
Variables: 

Pn 6 (0,Pn) - active generation at bus n 
5n £ ( — 00, 00) - reactive generation at bus n 
In t (0, 1.) - percentage of load served at bus n 
Variables from Model [2] (The Warm-Start LPAC Model) 
Maximize: 

In 



E 



(IV 4 



n£JV 

Subject to: 

Pn 



1) 



-Piln 
-Q l Jn - 



p g n VneJV 
-ql VneiV 



(IV 4 
(IV 4 



ql = Vn e N \ G 

n^m 

In = <j n m 
mS N 

Constraints from Model [2] (The Warm-Start 



-q A 



Vn e G 



(IV 4 



2) 
3) 

(IVftU) 



5) 

LPAC Model) 



maximum amount of load that can be dispatched. The model 
assumes that the loads can be shed continuously and that the 
active and reactive parts of the load should maintain the same 
power factor. The objective function (M4]l) maximizes the 
percentage of served load. Constraints and (Ivj4]3) set 

the active and reactive injection at bus n appropriately based 
on the decision variables for load shedding and generation 
dispatch. Constraint (!Vj4]4) ensures that reactive generation 
only occurs at generator buses and Constraint (]Vj4]5) now 
defines q n for generator buses as well. 

Since it reasons about reactive power and voltage magni- 
tudes, Model [4] can be further enhanced to impose bounds on 
these values. As we will show, such bounds are often critical 
to obtain high-quality solutions in power restoration contexts. 
If a reactive generation bound qs is supplied, this model can 
be extended by adding the constraint, 

q 9 n < & Vn € N. 

Voltage magnitude limits can also be incorporated. Given 
upper and lower voltage limits \V\ and \V\, the constraint 

\V\< l.O + n < \V\ VneN. 

may be used to enforce bounds on voltage magnitudes. The 
experimental results study the benefits of the LPAC model, 
suitably enhanced to capture these extensions, for power 
restoration. They compare a variety of linear models including 
the LDC model, the LPAC model, and enhancements of the 
LPAC model with additional constraints on reactive power and 
voltage magnitudes. 



Table VIII studies the applicability of various linear power 
models for network restoration on the IEEE30 benchmark. 
1000 line outage cases were randomly sampled from each 
of the TV- 3, N - 4, N - 5, . . . , TV - 20 contingencies. Each 
contingency is solved with a linear power model (e.g., the 
LDC model or the LPAC model), whose solution is used as 
a starting point for the AC model. The performance metric 
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TABLE VIII 

Power Restoration: Achieving AC Feasibility From Different 
Models. 



Scenario 


LDC 


LPAC 


LPAC-R 


LPAC-R-V 


N-3 


998 


999 


1000 


1000 


AT— 4 


999 


1000 


1000 


1000 


N—5 


987 


994 


1000 


1000 


N—6 


507 


594 


903 


1000 


N—7 


738 


856 


973 


974 


N—8 


949 


996 


1000 


1000 


N—9 


847 


932 


1000 


1000 


N—10 


219 


452 


992 


999 


N—ll 


726 


972 


1000 


997 


N—12 


491 


779 


998 


999 


Af-13 


444 


617 


983 


991 


Af-14 


545 


637 


998 


1000 


Af-15 


1000 


1000 


1000 


1000 


Af-16 


989 


1000 


1000 


1000 


Af-17 


1000 


1000 


1000 


1000 


Af-18 


969 


1000 


1000 


1000 


Af-19 


999 


1000 


1000 


1000 


Af-20 


1000 


1000 


1000 


1000 



TABLE IX 

Power Restoration: Average Load Shedding (% of Total Active 

Power). 



Scenario 


LDC 


LPAC 


LPAC-R 


LPAC-R-V 


Af-3 


3.23 


6.193 


15.94 


15.99 


Af-4 


2.827 


3.55 


9.183 


9.278 


Af-5 


0.9562 


2.204 


8.027 


8.082 


Af-6 


5.805 


6.149 


9.287 


9.921 


Af-7 


1.506 


5.709 


19.43 


19.46 


Af-8 


13.78 


18.54 


27.37 


27.39 


Af-9 


15.92 


24.27 


42.76 


42.78 


AT— 10 


29.23 


24.02 


32.09 


32.25 


AT-11 


25.18 


25.24 


42.62 


42.63 


AT- 12 


36.35 


28.25 


38.03 


38.82 


Af-13 


40.1 


32.56 


38.3 


38.66 


Af-14 


39.98 


36.9 


40.45 


40.67 


Af-15 


81.91 


81.92 


81.92 


81.92 


Af-16 


86.21 


86.31 


86.32 


86.32 


Af-17 


89.89 


89.89 


89.89 


89.89 


Af-18 


88.26 


88.3 


88.32 


88.32 


Af-19 


85.9 


86.13 


86.13 


86.13 


Af-20 


86.2 


86.37 


86.38 


86.38 



is the number of cases where the AC solver converges, as a 
good linear model should yield a feasible generation dispatch 
with a good starting point for the AC solver. To understand the 
importance of various network constraints, four linear models 
are studied: the traditional LDC model; the LPAC model; the 
LPAC model with constraints on reactive generation (LPAC- 
R); and the LPAC with constraints on reactive generation and 
voltage limits (LPAC-R-V). The number of solved models for 



each of the contingency classes is presented in Table VIII 
The results indicate that a traditional LDC model is overly 
optimistic and often produces power dispatches that do not 
lead to feasible AC power flows (the iV— 10 and N— 13 are 
particularly striking). However, each refinement of the LPAC 
model solves more contingencies. The LPAC-R-V model is 
very reliable and is able to produce feasible dispatches in 
all contingencies except 40 . This means that the LPAC-R- 
V model solves 99.76% of the 17,000 contingencies studied. 
Table |IX| depicts the load shed by the various models. For large 
contingencies, the LPAC-R-V model not only provides good 
starting points for an AC solver but its load shedding is only 



slightly larger than the (overly optimistic) LDC model. These 
results provide compelling evidence of the benefits of the 
LPAC model for applications dealing with situations outside 
the normal operating conditions. In addition, Model [4] can 
replace the LDC model in power restoration applications (e.g., 
P31l . [46]) that are using MIP models to minimize the size of 
a blackout over time. 

B. The Capacitor Placement Problem 

The Capacitor Placement Problem (CPP) is a well-studied 
application ll28l . Il29l . 11301 and many variants of the problem 
exist. This section uses a simple version of the problem to 
demonstrate how the LPAC model can be used as a building 
block inside a MIP solver for decision-support applications. 

Informally speaking, the CPP consists of placing capacitors 
throughout a power network to improve voltage stability. The 
version studied here aims at placing as few capacitors as 
possible throughout the network, while meeting a lower bound 
| V\ on the voltages and satisfying a capacitor injection limit q c 

and reactive generation limits q n (n g G). Model [^presents a 
CPP model based on the cold-start LPAC model. For each bus 
n, the additional decision variables are the amount of reactive 
support added by the capacitor q c n and a variable c n indicating 
whether a capacitor was used. 

The objective function (N^5]l) minimizes the number of 
capacitors. Constraints (1V(5]2) ensure the voltages do not 
drop below the desired limit and do not exceed the preferred 
operating condition of 1.05 Volts p.u. Constraints (1V(5]3) link 
the capacitor injection variables with the indicator variables, a 
standard technique in MIP models. Constraints (Ivj5]4) ensures 
each generator n e G does not exceed its reactive generation 
limit q^. and constraints (]vj5]5) defines the reactive power for 
generators. Lastly, Constraints (N^5]6) redefines the reactive 
power equation to inject the capacitor contribution q c . The 
remainder of the model is the same as Model [3] (the cold-start 
LPAC model). 

The CPP model was tested on a modified version of 
the IEEE57 benchmark. All of the IEEE benchmarks have 
sufficient reactive support in their normal state. To make an 
interesting capacitor placement problem, the transformer tap 
ratios are set to 1.0 and existing synchronous condensers are 
removed. This modified benchmark (IEEE57-C) has signif- 
icant voltage problems with several bus voltages dropping 
below 0.9. By design, a solution to Model [5] satisfies all of the 
desired constraints. However, Model [5] is based on the LPAC 
model and is only an approximation of the AC power flow. To 
understand the true value of Model [5] we solve the resulting 
solution network with an AC solver and measure how much 
the constraints are violated. Table |X] presents the results of 
Model [5] on _benchmark IEEE57-C with q c = 30 and various 
thresholds \V\. The table presents the following quantities: 

The minimum desired voltage |V|; The worst violation of 
the voltage lower-bound min(|T/|); The worst violation of the 
voltage upper bound max(|V|); The worst violation of reactive 
injection upper bound max(g„); The number of capacitors 
placed c n\ and the runtime of the MIP to prove the optimal 
placement solution. The table indicates that the CPP model is 
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Model 5 A MIP for the Capacitor Placement Problem 
Inputs: 

q n - injection bound for generator n 

q c - capacitor injection bound 

\V\ - minimum desired voltage magnitude 

Inputs from Model [5] (The Cold-Start LPAC Model) 
Variables: 

q n £ (0,9 C ) - capacitor reactive injection 
c n 6 {0, 1} - capacitor placement indicator 
Variables from Model [5] (The Cold-Start LPAC Model) 
Minimize: 

ngiV 

Subject to: 

\V\ < 1.0 + <f> n < 1-05 VneiV 

<?T< M c " 
Qn <Qn Vn e G 

n^m 

9™ = + Vn 6 G 

m£JV 

In — In = 'Y2 ?nm + <?rwn 
m£JV 

Constraints from Model [5] (The Cold-Start LPAC Model) except 



(M5jl) 



Vn£JV:n^sAn^G (TV 5 



(V5 2) 
(V5 3) 
(V@4) 

(M5 5) 



6) 
^5) 



TABLE X _ 
Capacitor Placement: Effects of \V\ on IEEE57-C, ~cf 



■ 30 MVar 



\v\ 


min(|V|) 


max(|V|) 


max(g n ) 




Time (sec.) 


0.8850 


0.000000 


0.0 


0.0 


1 


1 


0.9350 


0.000000 


0.0 


0.0 


3 


8 


0.9600 


0.000000 


0.0 


0.0 


5 


156 


0.9750 


0.000000 


0.0 


0.0 


6 


177 


0.9775 


0.000000 


0.0 


0.0 


6 


139 


0.9800 


0.000000 


0.0 


0.0 


6 


75 


0.9840 


-0.000802 


0.0 


0.0 


7 


340 



extremely accurate and only has minor constraint violations on 
the lower bounds of the voltage values. It is important to note 
that, although the CPP model can take as long as five minutes 
to prove optimalit}]^] it often finds the best solution value 
within a few seconds. The voltage lower bound approaches 
the value of 0.985, which is the lowest value of the voltage- 
controlled generators in the benchmark. These results remain 
consistent for other voltage bounds. 

Once again, the CPP model indicates the benefits of the 
LPAC approximation for decision-support applications that 
need to reason about reactive power and voltages. 

VI. Related Work 

Many linearizations of the AC power flow equations have 
been developed H, (6), Q, ED, EE E2), ED- Broadly, 
they can be grouped into iterative methods [6|, |16|, [59| and 
convex models Bl. 171. l35l, IT391. 

Iterative Methods: Iterative methods, such as the fast- 
decoupled load flow 1 59 1, significantly reduce the computation 
time of solving the AC equations and demonstrate sufficient 
accuracy. Their disadvantage however is that they cannot be 
efficiently integrated into traditional decision-support tools. 
Indeed, MIP solvers require purely declarative models to 
obtain lower bounds that are critical in reducing the size 
of the search space. Note however that, modulo the linear 

3 It is of course only optimal up the quality of the LPAC approximation. 



approximations, the LPAC model can be viewed as solving a 
decoupled load flow globally. The key differences are: 

1) Because the model forms one large linear system, all 
of the steps of the decoupled load flow are effectively 
solved simultaneously; 

2) Because the formulation is a linear program, the values 
of p and q can now be decision variables, and bounds 
may be placed on the line capacities, voltage magni- 
tudes, and phase angles; 

3) The model may be embedded in a MIP solver for making 
discrete decisions about the power system. 

The second and third points represent significant advantages 
over the fast-decoupled load flow and other iterative methods. 

Convex Models: Although many variants of the LDC 
model exist, few declarative models incorporate reactive flows 
in cold-start contexts. To our knowledge, three cold-start ap- 
proaches have been proposed: (1) a polynomial approximation 
scheme ll39l , (2) a semi-definite programming relaxation Q, 
and (3) a voltage-difference model l35l . 

The polynomial approximation has the advantage of solving 
a convex relaxation of the AC power equations but the number 
of variables and constraints needed to model the relaxation 
"grows rapidly" 11391 and only second-order terms were con- 
sidered. The accuracy of this approach for general power 
flows remains an open question: Reference 11391 focuses on 
a transmission planning application and does not quantify the 
accuracy of the approximation relative to an AC power flow. 

The semi-definite programming (SDP) relaxation J7J has 
the great advantage that it can solve the power flow equations 
precisely, without any approximation. In fact, reference Q 
demonstrated that the formulation finds the globally optimal 
value to the AC optimal power flow problem on a number of 
traditional benchmarks. However, recent work has shown this 
does not hold on some practical examples [60]. Computation- 
ally, SDP solvers are also less mature than LP solvers and their 
scalability remains an open question |61|. Solvers integrating 
discrete variables on top of SDP models [62 1 are very recent 
and do not have the scientific maturity of MIP solvers l48l . 

The voltage-difference model [35] has a resemblance to a 
model combining the equation 



v h =IK, 



h\2. 



\v*\\v*\9nm - \v£\\v£\b nm (e° n - e°j 



with Equation ( 17 1. However, it makes a fundamental assump- 
tion that all voltages are the same before computing the voltage 
differences. In practice, voltage-controlled generators violate 
this assumption. On traditional power system benchmarks, we 
observed that the voltage-difference formulation had similar 
accuracy to the LDC model. 

VII. Conclusion 

This paper presented linear programs to approximate the AC 
power flow equations. These linear programs, called the LPAC 
models, capture both the voltage phase angles and magnitudes, 
which are coupled through equations for active and reactive 
power. The models use a piecewise linear approximation of 
the cosine term in the power flow equations. The cold-start 
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Quality of a Piece-wise Linear Cosine Approximation 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Radians 



Fig. 4. A Piecewise-Linear Approximation of Cosine using 7 Inequalities. 

and warm-start models use a Taylor series for approximating 
the remaining nonlinear terms. 

The LPAC models have been evaluated experimentally over 
a number of standard benchmarks under normal operating con- 
ditions and under contingencies of various sizes. Experimental 
comparisons with AC solutions on a variety of standard IEEE 
and MatPower benchmarks shows that the LPAC models are 
highly accurate for active and reactive power, phase angles, 
and voltage magnitudes. The paper also presented two case 
studies in power restoration and capacitor placement to provide 
evidence that the cold-start and wam-start LPAC models can be 
efficiently used as a building block for optimization problems 
involving constraints on reactive power flow and voltage 
magnitudes. As a result, the LPAC models have the potential 
to broaden the success of the traditional LDC model into 
new application areas and to bring increased accuracy and 
reliability to current LDC applications. 

There are many opportunities for further study, including the 
application of the LPAC models to a number of application 
areas. From an analysis standpoint, it would be interesting 
to compare the LPAC models with an AC solver using a 
"distributed slack bus". Such an AC solver models the real 
power systems more accurately and provides a better basis for 
comparison, since the LPAC models are easily extended to 
flexible load and generation at all buses. 

Appendix A 

A Linear Programming Approximation of Cosine 

The convex approximation of the cosine function is imple- 
mented through a piecewise linear function that produces a 
linear program in the following way. The modeler selects 
a desired domain (I, hj^ and a number segments s. Then s 
tangent inequalities are placed on the cosine function within 
the provided domain to approximate the convex region. Figure 
|4] illustrates the approximation approach using seven linear 
inequalities. The dark black line shows the cosine function, 
the dashed lines are the linear inequality constraints, and 
the shaded area is the feasible region of the linear system 
formed by those constraints. The inequalities are obtained from 
tangents lines at various points on the function. Specifically, 

4 The domain should not exceed the range (—it/2, tt/2) to ensure convexity. 
In practice, £>° — 6^ is typically very small and a narrower domain is 
preferable. 



PWL<COS>(x £?s , x, I, h, s) 

1 postix^ > cos W-p(0 (x -l)+ cos(/)) 

2 inc<- (h- l)/(s + 1) 

3 a <— I + inc 

4 for i e l..s 

5 do f a <— cos(a) 

6 s a < sin(a) 

7 postix^ < s a x - s a a + f a ) 

8 a <— a + inc 

Fig. 5. Generating Evenly Spaced Piecewise-Linear Approximations. 

for an x-coordinate a, the tangent line is y = — sm(a)(x—a) + 
cos(a) and, within the domain of (— n/2, tt/2), the inequality 

y < — sin(a)(x — a) + cos(a) Va 

holds. Figure [5] gives an algorithm to generate s inequalities 
evenly spaced within (l,h). In the algorithm, a; is a decision 
variable used as an argument of the cosine function and x^Ss is 
a decision variable capturing the approximate value of cos(x). 
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